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Vibration-Translation Energy Transfer in Anharmonic Diatomic Molecules: 
II. The Vibrational Quantum-Number Dependence 


Robert L. McKenzie 

Ames Research Center, NASA, Moffett Field, California 94035 


I. INTRODUCTION 


Modern applications of vibrational energy transfer frequently deal 
with processes that depend on the details of excitation for specific vi- 
brational states. Infrared gas lasers are a primary example . 1 Of the 
several mechanisms Influencing the population of a particular vibrational 
state, the collisional conversion of vibrational energy to translation 
can be an important aspect. For example, the probability of vibration- 
translation (V-T) energy transfer is well-known to Increase as the quan- 
tum number of the initial state increases. Thus, even in situations 
where the V-T process may be Insignificant to the kinetics of lower vi- 
brational levels, it can dominate the flow of vibrational energy from 
upper levels and, in some cases, act as the primary path for vibrational 
energy loss from the system. The quantum-number dependence of V-T rates 
must therefore be considered in the analysis of most nonequilibrium pro- 
cesses where excess vibrational excitation has been produced. Unfortun- 
ately, very little quantitative information defining the V-T rate depen- 
dence on quantum number for even the simplest diatomic molecules is 
presently available from either experiment or theory. Experimental 
ground-state excitation rates have been obtained from measurements behind 
shock waves 2 or in fluorescence experiments^ for many years but the dif- 
ficulty of obtaining experimental V-T rates for molecules in well-defined 
excited vibrational states is Indicated by the sparsity of attempts. 
Numerous experimenters have recently measured the rates of vibration- 
vibration (V-V) energy exchange between pairs of oscillators in excited 
states 1 *^ because the fast V-V transfer can easily be made a dominant 
mechanism; but to date, only one comprehensive sec of upper level V-T 
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rate measurements has been reported . 4 Even then, while the experiment 
was cleverly designed and carefully analyzed, the conditions were complex 
and the measurements required substantial correction to compensate for 
extraneous modes of energy transfer. Theoretical studies addressed to 
the analysis of initially excited oscillators have been similarly sparse. 
The usual theoretical interest has centered on harmonic oscillators ini- 
tially in the ground state . 6 Simple analytic versions of these theories 
have frequently been applied in kinetic models to oscillators in excited 
states but their suitability in that application has not been validated. 

The purpose of this study is to extend the present analytical situa- 
tion by calculating the vibrational quantum-number dependence of V-T rate 
coefficients using a sufficiently complete collision model, which is not 
constrained to first-order approximations of the oscillator motion or its 
interactions. However, the approach to be taken is still limited by 
pragmatic considerations. Accurate V-T rate coefficient calculations by 
any theoretical model are obviated by uncertainties in the shape and mag- 
nitude of the interaction forces between colliding pairs for all but a 
few simple cases. Thus, we can only examine the qualitative features 
that are not masked by interaction potential uncertainties. Furthermore, 
even an extended collision model must retain some approximations, parti- 
cularly regarding, the collision geometry, if it is to remain computation- 
ally practical in the prediction of thermally averaged rate coefficients. 
Hence, a complete quantum-mechanical description of vibrationally inelas- 
tic encounters is avoided, although such descriptions have been formulated 
and solved with all degrees of completeness . 7 " 12 Instead, attention is 
confined here to a semiclassical treatment that accurately reproduces all 



of the main characteristics of vibrational energy transfer to initially 
excited oscillators but may be further reduced to yield closed-form 
analytic solutions. The analytic solutions are of primary interest here 
because of their practical Importance in the numerically cumbersome analy- 
sis of macroscopic nonequilibrium processes in which rate information for 
several modes of energy transfer must be economically provided . 1 The 
complete semiclasslcal model, requiring numerical solution, is applied 
both to an examination of the qualitative nature of upper state transi- 
tions and as a basis for evaluating the accuracy of the analytical solu- 
tions . 

In the sections to follow, the features of the collision model that 
appear most important to the dynamics of a vibrationally excited oscil- 
lator are first discussed; followed by a description of several approxi- 
mations, each of which retains one or more of the features considered. 
Approximate values of the interaction potential parameters and their range 
of uncertainty are then estimated by comparing the predicted ground-state 
rate coefficient with a comprehensive set of experimc ■ i, values. Colli- 
sions of CO with He are chosen as the example because of the abundant 
data available. The Implied potential parameters are then used to com- 
pare the numerical model with the experimental excited-state rate coeffi- 
cients and with the analytic predictions. Finally, the effects of 
multiple-quantum transitions from excited states on a vibrational relax- 
ation process are. considered both for molecules like CO where the effect 
is secondary and for molecules like the halogens where the effect can be 
dominant. 
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II. THE COLLISION MODEL 


A. Features influencing the excited state collision dynamics 

As the quantum number of the initial oscillator state is increased, 
several aspects influencing the oscillator dynamics and its interactions 
with the incident particle gain Increasing Importance. For example, the 
wavefunctions describing vibrationally excited eigenstates become more 
extended in the oscillator coordinate. Consequently, when the oscillator 
is distorted by a collision, the wavef unction overlap is greater not only 
with adjacent eigenstates but with more remote states as well. This 
feature is reflected by the increased magnitude of the matrix elements 
dynamically coupling the eigenstates which, in turn, accounts for the 
greater probability of V-T energy transfer through both single- and 
multiple-quantum transitions. Furthermore, the increased coupling of 
nonadjacent states during the collision can affect the final occupation 
of states adjacent to the initial state and thereby influence the rate 
of single-quantum transitions. Thus, a calculation of the oscillator 
dynamics from an excited initial state must include multiple-state 
interactions at collision energies where they are normally unimportant 
for oscillators in the ground state. 

The degree with which multiple-quantum transitions Influence the 
oscillator dynamics during a collision depend, in part, on the form of 
the interaction potential. A common practice, often used to simplify the 
analysis of ground-state oscillators, is to consider the oscillator mo- 
tion to be small compared to the range of interaction and linearize the 
interaction potential in the oscillator coordinate. In a harmonic 
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oscillator, this treatment has the effect of equalizing all of the diag- 
onal matrix elements and forbidding multiple-quantum transitions. The 
occupation of nonadjacent oscillator states is then possible only through 
a sequence of single-quantum steps during the collision. Nonlinear inter- 
action terms remove these restrictions and modify final state occupations 
in two related ways. First, ail of the nonadjacent states are directly 
coupled, thereby increasing their accessibility. Second, the diagonal 
matrix elements are no longer equal, leading to additional phase distor- 
tions in the quantum-mechanical oscillator motion that modify the proba- 
bility of transition. The additional phase shifts depend on the product 
of the difference between diagonal matrix elements and the strength of 
the interaction. They appear explicitly in a semiclassical impact param- 
eter treatment described by Bates 13 and applied to anharmonic oscillators 
by Mies. 11 * The formulation is reviewed in a subsequent section of this 
paper. 

All of the foregoing effects are amplified when oscillator anharraon- 
iclty is included. Nonadjacent states become coupled even for linearized 
interactions and the larger difference between the diagonal matrix ele- 
ments creates phase distortions that can become a significant fraction 
of the unperturbed oscillator period. Mies 11 * has shown the influence on 
transition probability predictions to be large even for oscillators ini- 
tially in the ground state. A second, and in some cases greater, effect 
of anharmonicity is its influence on the variation of eigenenergies with 
quantum number. Since transition probabilities and the related rate coef- 
ficients are well-known to depend on the amount of energy transferred, a 
lowest-order effect of oscillator anharmonicity may be demonstrated by 
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simply Inserting anharmonic oscillator eigenenorgies Into a harmonic os- 
cillator theory such as that given by Schwartz et aJL. 15 The results devi- 
ate substantially from the simple Landau-Teller relation for the rate 
coefficients given by 


fc . (T) 
m,m-L 


° mk l,0< T > 


( 1 ) 


where k m (T) denotes the rate coefficient for transitions from state 
m :o m-i and is a function of the kinetic temperature T. However, 
the simple aid hoc insertion of anharmonic eigenenergles into a harmonic 
escalator model is not always a sufficient means of accounting for an- 
harmonicity. The influence of anharmonicity on the interaction matrix 
elements, which in turn effects both the magnitude and phase of the oscil- 
lator motion, is often so great that an anharmonic oscillator model must 
be used from the start. Fortunately, oscillator anharmonicity and non- 
linear interaction potentials present only a slight increase in computa- 
tional difficulty, particularly if a Morse oscillator and an exponential 
form of the interaction are adopted. The necessary matrix elements are 
then conveniently expressed in closed algebraic form 9 > 1, *» 1 ^ just as they 
are for harmonic oscillators. 

Finally, an oscillator potential creating anharmonicity also admits 
to the existence of continuum states. We shall neglect their contribution 
to the energy transfer process, however, since they are energetically 
inaccessible by a large margin for the combinations of collision energies 
and initial states to be considered here. Although no evaluation has 
been made of their effects, the occupation of continuum states is presumed 
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to be os Bmall as the nearby bound states, and no bound states near the 
continuum were iound to influence the dynamics of any states at tne quan- 
tum levels of Interest. 

B. The semiclassical numerical model 

To obtain V-T rate coefficients, we calculate the associated transi- 
tion probabilities including oscillator anharmonlcity and nonlinear 
interactions by treating the oscillator quantum-mechanlcally but calcu- 
lating the collision trajectory classically. The trajectory is further 
constrained to collinear encounters. The unperturbed molecular wave- 
functions are those of a quantized Morse oscillator and the Interaction 
potential is composed of exponentially repulsive and, if desired, attrac- 
tive terms. Formulation of the numerical model has been described in 
detail previously 16 but, in brief, an arbitrary multitude of coupled 
Morse eigenstates are included in a time-dependent numerical solution of 
the Schrddlnger equation. Convergence is ensured by including a suffi- 
cient number of states above and below those of interest. The classical 
trajectory is coupled to the oscillator motion in a rigorous manner (with- 
in the semiclassical framework) guided by Ehrenfest's theorem. In effect, 
the classical path is influenced by the oscillator compression and recoil 
during the encounter. This feature successfully extends the semiclassi- 
cal method to encounters where the average size of the oscillator is 
severely perturbed during the collision (e.g. , when a light oscillator 
nucleus is struck by a heavy collision partner, as in H2-He or CO-Ar 
collisions). The predicted transition probabilities have been shown 16 
to reproduce the results from an equivalent collinear quantum-mechanical 
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calculation 9 over a wide range of conditions. We find, however, that al- 
though the semiclassical method Is conceptually simpler and more easily 
reduced to an analytic formula, numerical solutions of the complete semi- 
clasBicul model appear to be no more economical than exact quantum- 
mechanical solutions obtained using modern algorithms optimized for the 
roblem .®> 9 One of the penalties of using a semiclassical approximation 
is Chat total energy is not conserved, but the effects of that omission 
are easily and accurately compensated for by interpreting the relative 
collision energy or velocity as an average of the known initial and final 
values. A far mo.^e severe limitation of the semiclassical theory is its 
incomplete treatment of the interaction when the oscillator is very 
heteronuclear 1 6 (e.g., the hydrogen halides;. Such cases are avoided 
here and have presented numerical difficulty in exact treatments . 9 

The implications introduced by a restriction to collinear encounters 
are not as well understood but the restriction is necessa:y if the 
quantum-number dependence of thermally averaged rate coefficients is ever 
to be obtained in a reasonable computing time. Clearly, a more realistic 
approach would include a three-dimensional collision geometry in which 
simultaneous rotational transitions are coupled with the vibrational 
motion. Considerable activity along these lines is currently evi- 
dent 7-1 2 » 17,18 but the large number of rotational states that become 
accessible at collision energies sufficient to cause vibrational transi- 
tions would make our objective impractical for all but a few special 
molecules, like H2. On the other hand, as long as the rotational eigen- 
energies of the undisturbed molecule are well-described by a rigid-rotor 
model vsuggesting that the rotational and vibrational motion are separ- 
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able), the disparity between a collinear and a three-dimension! theory 
is not expected to be very sensitive to the Initial vibrotin ; al quantum 
number. By normalizing the collinear predictions accord! g to the ratio 
m _]/rok^ 0 anc * Qv °iding the prediction of absolute .ate coefficients, 
much of the absolute error associated with the co’ilnear restriction 
will hopefully be nullified. Such a ratio alpj absorbs the lowest-ordcr 
quantum-number dependence suggested by Eq, (1). 


C, Thermally averaged rate coefficients from a collinear scmiclassical 
model 

With the possible exception of molecular beam analyses, the appli- 
cations of an inelastic collision model usually require results in the 
form of a thermally averaged rate coefficient. A general formulation of 
the averaging integral is well-known but here the restriction to collin- 
ear trajectories and the use of a semiclassical approximation require 
some special consideration. In general, the rate coefficient for a ki- 
netic temperature T may be written in terms of the energy parameter 

c = E / kT and an energy-dependent cross-section a (E ) , where n 
mm m,n m 

denotes the final quantum state and E^ is the relative kinetic energy 
before a collision with an undisturbed oscillator in a pure eigenstate 
m. The rate coefficient is then 19 


-c 


k ( T ) = C 
mn 


a (E )e e m de , 
mn m m m 

* A 


( 2 ) 


where the average thermal speed is C = (8kT/mj)"^ and u is the reduced 
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collision mass. A further requirement for the collision model is that 
it conform to the detailed balance relations. Originating with tW 
reciprocity theorem, the requirements of detailed balance propagate 
through three levels of microscopic detail, giving the general physical 
relations for spinless nondegenerate collision partners as 


P 

mn 


< E m> 


P (E ) 
nm n 


(3a) 


E 0 (E ) « E a (E ) , 
m mn' nr n nm n 


-ftu> /kT -fun /kT 

k ( T ) e m = k (T)e 
mn nm 


(3b) 


(3c) 


where P is the transition probability from state m to n and tiw 
mn 111 

is the oscillator energy of state m. 

The collinear collision geometry produces semiclassical transition 
probabilities that behave according to Eq. (3a) but the restriction to a 
zero impact parameter leaves the cross-section required by Eq. (2) unde- 
fined. One common solution is to adopt an effective hard-sphere cross- 
section o and compute the Inelastic cross-section according to 


o (E ) = o P (E ) 
mn nr o mn m 


(4a) 


and 


a (E ) = a 'P (E ) . 
nm n o nm' n' 


(4b) 


Equation (3b) requires that 

°o’ = [1 + h( "n'“m )/E n lo o ’ 
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thus suggesting that the "hard-sphere” size must depend on the collision 
energy and transition in questionl This contradictory result is a conse- 
quence of the collineur approximation but the error is negligible when 

ki(w -to )|/E « I. In circumstances where the ratio approaches unity, 

1 n m 1 n 

the transition probability is typically so small that the integral in 
Eq. (2) is unaffected. 

Equation (2) must be further modified to compensate for the lack of 
energy conservation inherent in the semiclassical approximation. This 
discrepancy is easily and accurately corrected by interpreting the semi- 
classical relative collision energy E or speed u as an average of the 
initial and final values. Reference 16 demonstrates that while the cor- 
rection can be large, t‘he method of averaging has no apparent effect on 
the outcome for vlbratlonally Inelastic collisions at all energies from 
threshold up to the limits of practical interest. For convenience, we 
use an arithmetic energy average. Denoting the total energy as E^, the 
semiclassical approximation is brought into close agreement with an equiv- 
alent quantum-mechanical calculation by the interpretation 

E = E T - •n(ai jn + m n )/2 . (5) 

The combination of Eqs. (2), (4), and (5) then gives the thermal averag- 
ing prescription for a colllnear semiclassical collision model as 


k (T) = 
mn 


o c,>W 2kT 

o 


p (E) 

mn 


A Hnn 
. £ + 2kT J 


e e d c 


( 6 ) 


where e - E/kT and to = w - to . To make the satisfaction of Eq. 

mn m n 

(3c) by Eq. (6) more obvious, the lower Integration limit in Eq. (6) has 


12 



been set to zero even though the independent variable transformation 

from e to e via Eq. (5) produces a limit of ± l-fiu) /2kT | , depending 
m mn 

on the sign of w , The negative limit may clearly be reset to zero 
mn 

but, even when the limit is positive, the probability threshold is near- 
ly twice the limit, so that again setting it to zero has no effect on 
the integral. 

HI. ANALYTIC APPROXIMATIONS 

Of the many analytic approaches appearing in the literature (see 
Ref. 6 for a partial summary), three that stand out In their application 
for estimating the V-T rate coefficient variations with quantum number 
are (a) the semiempirical formulas for Morse oscillators of Keck and 
Carrier, 21 (b) the perturbation treatment of Morse oscillators developed 
by Mies, 11 * and (c) the exact solution to a linearly forced harmonic 
oscillator obtained by Kerner. 22 Each approach retains one or more of 
the aspects of special interest to this application. They share the 
common feature that all incorporate collinear collision geometry and 
are all based on an exponentially repulsive interaction potential (later 
referred to as potential I) of the form 

V I (y) = Ae" y/L , (7) 

where the coordinates are defined in Fig. 1 and transformed according to 
y = jc - yr. The mass ratio y = m^/tm^ + m^) locates the mass center of 
the molecule and L is an adjustable range parameter. 
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The Keck-Carrier formula for anharmorlc oscillators 


The formula obtained by Keck and Carrier 21 comes from an adaptation 
of the distorted-wave harmonic oscillator theoi/ of Schwartz et_ al. 15 for 
a Morse oscillator. It includes an empirical fit to the numerical solu- 
tion of an integral equation for the "adiabaticity factor" and provides 
a particularly simple formula for estimating single-quantum transition 
rates from an arbitrary initial state. Keck and Carrier 21 made no claim 
for the suitability of their formula in applications beyond a demonstra- 
tion of the role of vibrational nonequilibrium in a dissociating gas; 
but the formula was subsequently applied by Bray 23 in a pioneering and 
detailed calculation of a vibrational relaxation process for anharmonic 
oscillators, apparently because of its simplicity and for lack of a bet- 
ter estimate. For similar reasons, the Keck-Carrier formula has since 
gained widespread use in the detailed analysis of upper state kinetics 
in lasers. 1 Its consideration here is motivated primarily by the number 
of kinetic models that incorporate it. The Keck-Carrier formula can be 
written in a form similar to Eq. (1) as 23 


1-Xj F, 


k , (T) - m ~ k. ,, (T) , 

m,m-l 1 - mx o Fj 1,0 


where F is obtained from the empirical formula 


F = if 3 - a W3 | e W3 (9) 

in which 

n ** - w .LCii/akT)* 1 . UO) 

m,m-i 
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The transition frequency u . = tu - to . is computed for a Morse 

m,m-l m m-i 

oscillator from 

to *» it) [(m + *s) - x (w + H) 2 ] (11) 

me e 

where id is the fundamental oscillator frequency and x_ is the an ~ 

G ° 

harmonic correction. 


The Mies perturbation solution for anharmonic oscillators 


The closest approximation to the numerical model used here is a 

14b 

semiclassical first-order perturbation treatment developed by Mies. 

It properly includes the effects of anharmonicity but, by the nature of 
first-order methods, It neglects the influence of states other than the 
designated initial and final states. Furthermore, to obtain an analyti- 
cal solution the classical path must be computed independently from the 
motion of the oscillator. The theory is therefore applicable only to 


single-quantum transitions in which the transition probabilities are 
small compared to unity. The independent classical path further restricts 
its application to nearly homonuclear oscillators such as CO (and of 
course, all homonuclear molecules) colliding with atomic particles of 
lighter mass than either of the molecular nuclei. 16 The appearance of 
a probability maximum signals the failure of the theory. 16 In spite of 
these shortcomings, we shall see that Mies' solution still provides a 
more useful approximation of the numerical predictions than the other 
analytic formulas investigated. A convenient form of Mies' result for 


the transition probability from state m to n, where n = m± 1, can be 
14b, 16 

written 


15 


12 


P 

mn 


(E) 


E3 


^mn 2 ttruLu 

V fi sinh(itg) 
mm 


‘K-g,*) . 


( 12 ) 


where E is the relative collision energy in a mass-centered reference 
_ — U 

frame and u = (2E/p) is the corresponding speed. The other param- 
eters are g « Lw /ir and X ~ pLTT(\r - V m ) /-ftV . The series function 
mn nn mm mm 

^("g^) and the matrix elements V are defined in Appendix A. Note 

that equal diagonal matrix elements V = V lead to 4*(-£.0) = 1. 

mm nn ** ' * 

reducing Eq. (7) to the equivalent formula for a harmonic oscillator 

14b 

with a linearized Interaction. 

As with the numerical model, Eq. (12) produces energy-dependent 
transition probabilities while a temperature-dependent rate coefficient 


is desired. No analytic solution of the Integral Eq. (6) with P (E) 

mn 

given by Eq. (12) is apparent but a reasonably accurate technique 
(labeled "The Method of Steepest Descent") for obtaining an analytic 
approximation 6 is based on the well-defined maximum contained in the 
integrand of Eq. (6). The value of e at which the maximum occurs is 
determined primarily by the exponential arguments. 6 The remaining func- 
tion is slowly varying over the range of the integrand and may be evalu- 
ated at the single value e locating the peak. The exponential argu- 

P 

ment is then expanded i.d second order about the peak and the term inte- 
grated analytically. In this application, the notation is simplified 

U 

with the substitutions e = fill) /2kT and n = -m L(p/2kT) . The 

mn mn mn 

exponential nature of Eq. (12) is also simplified by noting that in the 
energy range where the perturbation analysis is applicable, the transi- 
tion period t = 2ir/m is typically less than the effective collision 
P mn 

period t £ = 2L/TT. Thus, 7rg = t c /t p >l and sinh(irg) » Equations 


1 6 


(6) and (12) are then combined to give 



V 2 1 

2 e 

k (T) ■ 0 c 


mn 

e 

mn o 

V € 

mm mnj 





(e+ l € mn l) e "fi+2i"le d - f (13) 


The integrand peak is located at 


e p = [y (™ mn L) 2 /2kT3 ^ . 


(14) 


Using the procedure described, the approximate solution to Eq. (13) be- 
comes; 


k (T) = 16 (3ir 3 )~^a C 


mn 


V e 


3) 


5E _E_ 

.V e 
(_ mm mn J 


c ^(e + |e | ) <t> 2 (-g ,X ) 
p v p 1 mn 1 °p’ p' 


-3e +e f 
*e P mn { 1 + erf 


OCpM )* 4 


(15) 


where g - e /tt and X = e 2 (V -V )/rre V . The error function 
“p p p p mm nn' mn mm 

in Eq. (15) is close to unity for most cases. Equation (15) has log k 
_U 

a T 3 as expected and satisfies Eq. (3c). The temperature at which a 
given collision speed is coincident with the peak of the integrand in 
Eq. (13) defines the most effective speed at that temperature; this tem- 
perature will also be useful and can be identified from Eq. (14) as 

T = yu 3 /2irkU jb . (16) 

p ‘ mn 1 

Comparisons of the approximate integration of Eq. (13) with exact numeri- 
cal integrations show that the approximate method is most accurate at low 
temperatures. The first-order perturbation formula, Eq. (12), is most 
accurate at low energies, thus further contributing to the accuracy of 
Eq. (15) at low temperatures. 


17 



c. 


Th e Kerner solution for linearly forced harmonic oscillators 


The final analytic formula to be considered is an exact solution 
obtained by Kerner 22 for a harmonic oscillator that undergoes a forcing 
function linear in the oscillator coordinate. That condition may be 
satisfied in situations where r/L « 1 in Eq, (7), The potential may 
then be linearized according to 


V(x,r) = Ae *^ L (1 - yt/h) . 


(17) 


Kerner* s solution was applied by Treanor 21 * in a semiclassical collinear 
approximation using Eq. (17). Within the framework of the collision 
model, the resulting formula exactly calculates the probability of transi- 
tions between arbitrary states with the interaction of all states included. 
Thus, it can be applied at high collision energies where the interactions 
of more than two states influence the oscillator dynamics. In spite of 
the approximate nature of the harmonic oscillator model, wherein direct 
multiple-quantum transitions and the unbalanced coupling of higher and 
lower states caused by anharmonicity are excluded, the Kerner solution 
remains useful because it offers the only analytic means of estimating 
transition probabilities at high energies. Examples will be shown where 
multiple-quantum transitions and oscillator anharmonicity are not domin- 
ant, allowing accurate prediction by the Kerner solution. 

Kerner 22 and Treanor 24 write the probability for transitions between 
two arbitrary states m and n as 


_ . , , -E_ „m+n 

P (E) = m! n! e ° E 
mn o 


J r 

l 

U~o L 


<-E 0 )* 


(m-j)! j! 


(n-j)l 


* 

r 


( 18 ) 
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where J is the lesser of the quantum numbers in and n. The parameter 
E q Is the energy absorbed by a classical harmonic oscillator divided by 
one quantum of vibrational energy. For a collinear collision and the 
interaction of Eq, (17), Rapp* 5 obtains 

E q = 2(2ho)Lyw) 2 e 2 ttwL/u^ u ^ . (19) 

In Eq. (19), is the reduced mass of the oscillator and n) is the 
oscillator frequency. The accuracy of the model, when applied to highly 
excited oscillators, is substantially improved if the effective oscilla- 
tor frequency is corrected for anharmonicivy for each initial state m 
according to m = 10^(1 - 2x g m) . Without the correction, the excited-state 
rate coefficients would simply behave according to the Landau-Teller 
relation, Eq. (1), at low temperatures where the effective values of E q 
are all less than unity and give k^ m _^/mk^ q< 1 for large E q . An 
inconvenience of the Kerner formula Is its incompatibility with the ap- 
proximate integration method of Eq. (6) for obtaining a rate coefficient. 
A simplified version of Eq. (18), assuming E q « 1 , permits an approxi- 
mate analytical solution. However, the calculations are then restricted 
to a thermal range where multiple-quantum effects are insignificant and 
the theory loses its advantages over perturbation solutions. In the com- 
parisons to follow, we have therefore resorted to a numerical integration 
of Eq. (6) when the Kerner solution is applied. 
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IV. COMPARISONS WITH CO-He EXPERIMENTS 


In this section, the ability of the theoretical model to reproduce 
experimental rate coefficients is tested. Unlike past comparisons of 
vibrational rate coefficients with theory, we now have access to at 
least one set of experimental values for excited initial states. 1 * 

To test the consistency of the theory and experiment for all vibra- 
tional states, however, the effective interaction range L and the 
hard-sphere cross-section o q are determined from the abundant collec- 
tion of measurements dominated by transitions between the ground state 
and first vibrational state. The interaction parameters required to 
match the ground-state experiments are then applied in comparisons with 
the excited-state rate measurements. 

A. Effective interaction parameters 

The computational convenience gained from the simplified interaction 
potential I, Eq. (7), justifies its use, but as a consequence of its sim- 
plified form, the predicted rate coefficients cannot be expected to repro- 
duce the experiments at all kinetic temperatures. Transitions induced 
in an oscillator depend to a large extent on the potential gradient near 
the distance of closest approach; while in a collinear collision, the 
distance of closest approach is determined mainly by the coordinate where 
the potential magnitude equals the initial kinetic energy of the colli- 
sion. The magnitude of a purely repulsive potential, such as Eq. (7), 
and that of a more realistic potential with an attractive well may be 
the same at the closest approach distance but have a significantly dit- 
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ferent gradient. Consequently, where collisions are averaged over a 
range of energies, the predicted variation of rate coefficients with 
kinetic temperature will be different for the two potentials. By match- 
ing theory and experiment in several thermal ranges, and by using more 
than one potential form, an indication of the degree of uncertainty in 
rate coefficients attributable to potential errors can be obtained. For 
that purpose, we consider a second potential given by 

VjjCy) = D e (y °" y ^ L - 2D e ^ y ° _y ^ 2L . (20) 

Potential II is a Morse-type interaction with an attractive well of depth 
~D at coordinate y . As with Eq. (7), the exponential form allows ma- 
trix elements to be calculated analytically. 

Predictions by the numerical enharmonic oscillator model with the 
oscillator initially in the first eigenstate m = 1 are compared with 
experiment in Fig. 2. When potential I, Eq. (7), is used, the rate 
coefficients are independent of the magnitude A, so that only the range 
L requires specification. Similarly, the predictions using potential 
II are independent of y Q but require both L and D to be specified. 

The value D/k = 100°K is representative of well depths inferred from 
viscosity measurements. 26 The two potential gradients are different by 
about 20% at closest approach for the typical conditions considered. 
Figure 2 demonstrates the expected results. No unique set of potential 
parameters reproduces the experiments over the complete thermal range 
but the more realistic potential II comes the closest. The required 
values of L fall between 0.02 nrn and 0.03 nm, depending on the thermal 
range considered. 
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As an interesting aside, note that the low temperature dv;parture of 

-U 

the experimental rate from a variation proportional to T 3 in also 
followed by the theory using simple repulsive potentials. As Shin 27 
points out, these low temperature departures do not necessarily depend 
on weak attractive forces normally omitted from the interaction poten- 
tial; they even occur with a repulsive potential when the thermal aver- 
aging integration is done accurately for low collision energies. We 
know, however, that real interaction potentials usually contain an 
attractive component and it will augment this low temperature behavior. 

B. Comparisons with excited-state rote measurements 

Normalized rate coefficients, predicted for initially excited CO at 

T = 300°K, are compared in Fig. 3 with the room temperature measurements 

of Hancock and Smith. 1 * The parameter k ,/mk. „ is much less sensi- 

m,m-i 1,0 

tive to interaction uncertainties than the absolute rate coefficients 
and varies in a simple, nearly linear manner with initial-state quantum 
number m. The nearly linear quantum-number dependence, increasing with 
m at room temperature, is predicted for all of the interaction poten- 
tials examined and is believed to be an accurate description of the real 
behavior. As Fig. 3 shows, the experimental excited-state values com- 
pare favorably in magnitude with the predictions, but their trend is 
inconsistent with a linear extrapolation to m e 1. A highly nonlinear 
extrapolation is contrary to any prediction of the collision model at 
any temperature. Although the collision model contains many simplifi- 
cations awaiting refinement, the behavior implied by the experimental 
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rates appears also to require further verification and extension. In 
the interim, the theoretical predictions of excited state rates seem to 
be qualitatively reasonable and self-consistent in spite of their quan- 
titative uncertainty. Unfortunately, their verification by experiment 
remains Inconclusive. 

V. AN EVALUATION OF THE ANALYTIC APPROXIMATIONS 

The computational expense of the numerical model makes it impracti- 
cal as a general means of estimating excited-state rate coefficients. 
Instead, it is used in this section as a basis for evaluating the more 
convenient but less complete analytic formulas. The predicted rate coef- 
ficient variations with quantum number for several models are illustrated 
in Fig. 4 for two extreme temperatures. Tne differences in the various 
models depend strongly on the kinetic temperature but they all predict 
a simple monotonic change with quantum number. The analytic approxima- 
tions are therefore more clearly evaluated by choosing the highest ini- 
tial quantum number of practical interest and then comparing the predic- 
tions for a range of temperatures. In the case of CO, Lord! et al. 1 ^ 
have shown that energy transfer from vibrational levels as high as the 

twentieth can influence the net energy balance in an electrically excited 

CO laser system. Choosing m = 20, the single-quantum rate coefficients 

predicted by all of the collision models are compared in Fig. 5. The 

independent parameter (hw e /kT)^ was chosen so that predictions by the 

Keck formula, Eq, (8), appear as a nearly straight line. A comparison 

of the rates from the numerical model using potentials I and II shows 

the moderate sensitivity of k . /mk. n to the form of the potential 

m ) m~ l 1 j U 
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for one potential range L at all temperatures. Not shown 1 b the great 
sensitivity of the magnitude of k^ j^/mk^ q t0 other potential ranges 
at any temperature. Note, however, that the qualitative nature of the 
predictions are undisturbed by the form of the potential and therefore 
considered realistic. As expected, the MieB solution, Eq. (15), accu- 
rately reproduces the numerical results at low temperatures but fails at 
higher temperatures where multiple-state interactions begin to affect 
the single-quantum transitions. The departure is signaled wi.en transi- 
tion probabilities approaching unity influence the thermal averaging 
integral, Eq. (6). Since CO is not very anharmonic, the Kerner harmonic 
oscillator model, Eq. (18), frequency-corrected for anharmonicity at 
m = 20, works well over the entire thermal range. Note that the anhar- 
monic correction must be included, however, as all predictions are sig- 
nificantly above the result stated by Eq. (1) for a single-frequency 
harmonic oscillator. Finally, Fig. 5 shows that the Keck formula, 

Eq. (8), is too crude on approximation for large initial quantum numbers. 

The degree of oscillator distortion caused by the collision of a 
light helium atom with a CO molecule has an insignificant effect on the 
classical trajectory. This fact is made evident by the small difference 
at low temperatures between the numerical model where the effect is in- 
cluded and the Mies solution where it is neglected. An example in 
which the coupling is larger is illustrated in Fig. 6 for C0(m=20)-Ar 
collisions. In this situation, none of the analytic models do well at 
low temperatures because the effects of oscillator distortion on the 
classical path modifies the transition probabilities even near threshold. 
The small corrections are then greatly amplified by the thermal averaging 
integral at low temperatures. 
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The Bmall anharmonicity of CO(x e a 0.0062) has influenced the pre- 
ceding examples mainly by altering the energy spacing between excited 
eigenstates. Anharmonicity also modifies the absolute magnitude of the 
rate coefficients but that effect is not apparent in the ratio 

m-l^ m * c l 0' e3cnin P^ e *- n which the anharmonicity is large is illus- 
trated in Fig. 7 for U 2 (m » 10)-He(x e ° 0.0268). in this case, the 
frequency-corrected harmonic oscillator model is Inaccurate at all tem- 
peratures. The large spacing between elgenenergies in U 2 suppresses the 
onset of multistate interactions at high temperatures, making the Mies 
solution an accurate reproduction of the numerical results over the 
entire thermal rangr, The difference in mass between the He and H nuclei 
produces only moderate coupling between the compressed oscillator and 
the classical path. 1 6 

As the preceding comparisons Indicate, one cannot generally choose 
a single analytic model for estimating excited-state rate coeffipients 
that is applicable to all collision pairs. The situations where a model 
should not be used are easier to identify. Clearly, the Keck formula, 

Eq. (8), is too approximate in all of the examples. The Kerner harmonic 
oscillator solution, Eq. (18), with anharmonlcity-corrected frequencies 
is reasonably accurate unless the anharmonicity is large. The Mies an- 
harraonic oscillator solution, Eq. (15), is a poor approximation when 
multiple-state interactions become important. Finally, no analytic model 
based on the semiclasslcal approximation will be realistic when the oscil- 
lator dynamics have a significant influence on the classical path of the 
incident particle. This restriction limits all of the me !ele considered 
to collision pairs in which the mass of the incident particle is not 
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significantly greater than the mass of the impacted nucleus and to oscil- 
lators that ..re not extremely heteronuclear . 

VI. MULTIPLE-QUANTUM TRANSITIONS 

In the preceding section, only transitions to an adjacent state have 
been examined. Here, we investigate the relative importance of multiple- 
quantum transitions, particularly for oscillators in highly excited 
states. The probabilities of multiple-quantum transitions are compared 
in Fig. 8 both for CO(m)-He collisions in which the oscillator is ini- 
tially in an excited state and in states near the ground state. The 
collision speeds contributing most to the thermally averaged rate coeffi- 
cient at a selected temperature are indicated by the effective temperature 
T . In the thermal range considered, multiple-quantum transitions to the 

p 

ground-state are always improbable compared to single-quantum transitions 
from the first vibrational level, but the situation is clearly different 
when the oscillator is initially in the twentieth quantum state. However, 
thermally averaging the transition probabilities in Fig. 8 reduces the 
apparent importance of multiple-quantum transitions in a relaxation pro- 
cess. Figure 9 illustrates the resulting rate coefficients for two poten- 
tial ranges, using potential I and values of o^ obtained from the experi- 
mental match in Fig. 2 at T = 1000°K. The amplified uncertainty caused 
by the Interaction potential and its influence on the implied value of o q 
is most obvious but the qualitative features are again consistent for both 
potential ranges. For oscillators like CO, multiple-quantum transitions 
provide a significant path for energy transfer only at very high tempera- 
tures, according to these predictions. 
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A temperature marking the onset of competitive multiple-quantum 
transitions Is the characteristic vibrational temperature of the oscil- 
lator, here defined as 0^ = -fiu^/k (for CO, 0 v = 3122°K.). An oscillator 
in which multiple-quantum transitions will dominate the relaxation pro- 
cess can then be identified if 6^ is small. One extreme example is 
Br£ for which 0^ = 465° K. Since the anharmonicity is also small in 
B r 2 (x “ 0*0033) , the Kerner harmonic oscillator model has been used to 

j obtain the Br 2 ~He rate coefficients displayed in Fig. 10. Two- and 

f 

three-quantum transitions from the tenth vibrational level are shown to 
be significant even at room temperature and the temperature dependence 
of the single-quantum rate (m-t n = 10-t 9) is Inverted by multiple-state 
interactions. The high probability of multiple-quantum transitions in 
this case contributes to the extremely fast and thermally insensitive 
relaxation rates measured in the halogens and destroys the concept of a 
single "relaxation time" 19 that is Independent of the nonequilibrium 
state of the process for molecules of this type. 


VII. CONCLUDING REMARKS 

We have relied on a collinear semlclassical model for vibrationally 
inelastic collisions entirely for pragmatic reasons. The collinear geom- 
etry affords an economically reasonable means of estimating V-T rate 
coefficients for excited molecules and the semiclassical approximation 
is easily reduced to practical analytic solutions. While these simpli- 
fications clearly obviate the quantitative accuracy of the calculations, 
n' serious omission is apparent that would modify their qualitative 
nature, even in the presence of uncertain interaction potentials. 
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Unfortunately, an attempt to confirm the predicted features through ex- 
perimental comparison was inconclusive. However, the experimental con- 
ditions that would test the model most severely can at least be identi- 
fied. For example, the choice of collision partner has a large influ- 
ence on the rate coefficient sensitivity to initial quantum number. 

Note the large deviations of k ,/mk, n from unity in Fig. 6 for 

m , m - 1 1 , U 

CO-Ar compared to those for CO-He in Fig. 5. Furthermore, the increased 
oscillator distortion caused by heavy atom impact requires a more complete 
description of the Interaction than needed for light atoms. From another 
viewpoint, the lesser sensitivity of some features of the prediction to 
uncertainties can guide the choice of experimental variables to be empha- 
sized. In particular, an apparently universal feature of the V-T 
excited-state rate predictions is their monotonic low-order variation 
with quantum number. Once this feature is confirmed, the experimental 
emphasis can be shifted to the less predictable variations with kinetic 
temperature. Finally, a comparison of the estimates using various poten- 
tial parameters suggest that a self-consistent set of experimental rates 
for both high and low initial quantum numbers contains much more informa- 
tion defining the interaction potential than ground-state rates alone. 

Comparisons of the analytic and numerical rate coefficients graph- 
ically delineate the suitable range of application for each analytic 
model. However, the utility of an analytic approximation can also depend 
on the physical properties of the application. For example, the Kerner- 
Treanor harmonic oscillator model, with anharmonically corrected frequen- 
cies, predicts the ratio k ,/mk. with surprising accuracy for 
r m,m-i l,u 

many molecules; but before the model can be economically applied, an 
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analytical solution to the thermal averaging Integral, Including Kerner's 
transition probability formula, awaits development. Even with that solu- 
tion in hand, one must be concerned with the effect of anharmonicity for 
each molecule treated by the model. On the other hand, Mies' solution 
for anharmonlc oscillators, Eq. (15), fails at high temperature. At 
those conditions, however, many nonequilibrium processes are insensitive 
to the V-T rates of excited states, either because the vibrational state 
population distribution is nearly Boltzmann or because the process is 
controlled by some separate energy transfer mechanism. At lower temper- 
atures, the model accurately deals with a broader range of oscillators 
because anharmonicity is rigorously included. Collision partners for 
which the theory fails are poorly treated by all the analytic solutions 
based on a seraiclassical approximation. Similarly, the frequently employed 
formula developed by Keck^ is useful because of its simplicity but the 
additional computation required by the Mies solution is not prohibitive. 

The series function •K-g.A) converges rapidly and the matrix elements 
may be computed in advance. 

The calculations of multiple-quantum transition rates from excited 
states validate the assumption most often made in kinetic models of non- 
equilibrium processes: they can usually be neglected. As before, at very 
high temperatures where multiple-quantum transitions become competitive, 
a nonequilibrium process is usually not controlled by excited-state V-T 
rates while ground-state transitions are still dominated by single-quantum 
steps. Molecules with closely spaced vibrational energy levels, such as 
the halogens, are notable exceptions requiring a more careful analysis. 



APPENDIX A 


The series function <1> (— g , X ) and the matrix element V for an 

mn 

anharmonic oscillator are given here in algebraic terms for use in Eqs. 
(12) and (15). A more complete description of their origin is given in 
Ref. 16. The series function is computed according to 


^ ° l V *" 1 > 

2,-1 * 


(Al) 


where k\ = 1, A 2 “ -g. The remaining coefficients are obtained by the 
recurrence formula 


A (2,- 1)A^ - “ZgA^ - a £ _ 2 


(A2) 


The matrix elements, defined by V - ^n|e |m^ , are most easily com- 

1. 

puted in terms of the parameters a = y/aL, 3 = x o -1 > and a = 2irm e (2|J o x e ) ‘ 


where all variables are defined in the main text. The result is 9 > ll+ > 16 


„ _ „a Vn T(B-n) r (-l)^ -1 " r(l + a + n-2) r(B-a-l-n) 

mn ~ P a m! U (n - ^) ! T(1 + a + n - m - £) V (3 - 2n + 2.) 


(A3) 


In Eq. (A3), r (z) is the gamma function with argument z and the normal- 
ization constants are 


N j = [a (3 - 1 - 2j)/r(B- j)]* 2 . 


(A4) 
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FIGURE CAPTIONS 


FIG. 1. 
FIG. 2. 


FIG. 3. 


FIG. 4. 


Collinear collision geometry. 

A comparison of experimental rate coefficients for CO(m = 1)-He 
transitions to the ground vibrational state with predictions 
from the numerical model 16 for anharmonic oscillators. The 
solid and long-short dash lines were computed using the repul- 
sive interaction potential I, Eq. (7). The short dash line 
was computed using the Morse interaction potential II, Eq. (20). 
Hard-sphere collision cross-sections were chosen for each po- 
tential to match the experiment at T = 1000°K. Experimental 
values are from: o Ref. 2, • Ref. 3a, A Ref. 3b, ^ Ref. 3c. 

A comparison of experimental rate coefficients at T ~ 300°K 
for C0(m)-He transitions from vibrational states m to m - 1 
with predictions from the numerical model using repulsive 
potential I, Eq. (7). The excited-state data are from Ref. A 
and have been normalized using the experimental q value 
of Milliken 2 » 3 (Fig. 2). 

The CO(m)-He rate coefficient dependence on quantum number 
predicted by several collision models. The solid lines repre- 
sent the anharmonic numerical model, 16 the long-short dash 
lines represent the Kerner harmonic oscillator solution, 22 
Eq. (18), and the dashed lines are from the formula of Keck, 21 
Eq. (8). The potential range L = 0.02 nm, was used in all 
cases. 
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FIG. 


FIG. 


FIG. 


FIG. 


FIG. 


FIG. 


5. A comparison of excited-state rate coefficients for C0(m°20)-He 
predicted by several collision models. The potential range was 
L = 0,02 ran i- 1 all cases. 

6. A comparison of excited-state rate coefficients for C0(m ® 20)-Ar. 
Potential I, Eq. (7), was used with L = 0.02 run in all coses. 

7. A comparison of excited-state rate coefficients for II? (m" 10)-He. 
Potential I, Eq. (7), was used with L » 0.02 nm in all coses. 

8. Multiple-quantum transition probabilities for C0(m)-He colli- 
sions using the enharmonic numerical model with potential I 
and L = 0.02 nm. The effective temperature T^ locates the 
mo3t effective collision speed contributing to the thermally 
averaged rate coefficient at the temperature designated. 

9. Multiple-quantum rate coefficients for C0(m)-He. Potential I 
was used in the enharmonic numerical model. The hard-sphere 
cross-section values o q for each potential range are those 
required to match the experimental rates in Fig. 2 at 

T - L000°K. 

10. Multiple-quantum rate coefficients for Br 2 (m)-He predicted 
using the Kerner harmonic oscillator solution, Eq. (18), with 
L = 0.02 nm. 
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FIG. 1. 


Col linear collision geometry. 
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FIG. 2. A .-ora pari sun of fxpiriwntal rate coefficients for C0(m»l)-He 
transitions to the ground vibrational state with predictions 
from the numerical model 1 1 for anharmonic oscillators. The 
solid and long-short dash lines were computed using the repul- 
sive interaction potential 1. Eq. (7). The short dash line was 
computed using the Norse interaction potential 11, Eq. (20). 
Herd-sphere collision cross-sect ions were chosen for each po- 
tential to match the experiment at T - 1000°K. Experimental 
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FIG. 3. A comparison of experimental rate coefficients ai T * 300°K 
for C0(m)-He transitions from vibrational states m to ■ - 1 
with predictions from the numerical model usinf repulsive 
potential I, Eq , (7). The excited-state data are from Ref. 4 
and have been normalized using the experimental k. value 
of Mil liken (Fig* 2)« * 
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FIG. A. The C0(m)-He rate coefficient dependence on quantum number 
predicted by several collision models. The solid lines 
represent the anharmonic numerical model, 1 * the long-snort 
dash lines represent the Ker*er harmoni- oscillator solution,'- 
Eq. (18). and the dashed lines are from the formula of Keck. ?I 
Eq. (8). The potential range L - 0.02 nm, was used in all 
cases. 
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FIG. 5. 


A comparison of excited-state rate coefficients for C0(m-20)-He 
predicted by several collision models. The potential range was 
L - 0.02 nm «n all cases. 
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FIG. 6. 


A comparison of excited-state rate coefficients for C0(m-20)-Ar 
Potential I. F.q. (7), was used with L * 0.02 nm in all cases. 
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FIG. 7. 


A comparison of e.-:cited-state rate coefficients for H 2 (m - i0)-He. 
Potential 1. Eq. (7), was used with L * 0.02 nm in all cases. 


TRANSITION PROBABILITY 


EFFECTIVE TEMPERATURE, T p , °K FROM EQ. ( 16 ) 



RELATIVE COLLISION SPEED, U, km/sec 

FIG. 8. Mu It iple-quaritum transition probabilities for CO(ra) - He colli- 
sions using the anharmonic numerical model with potential I 
and L « 0.02 run. The effective temperature T p locates the 
most effective collision speed contributing to the thermally 
averaged rate coefficient at the temperature designated. 
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Multiple-quantum rate coefficients for Br : (m)-He predicted 
using the Kerner harmonic oscillator solution, Eq. wi 
L * 0.02 nm. 


